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Abstract. - Stochastic differential equations in Hilbert space as random nonlinear modified 
Schrodinger equations have achieved great attention in recent years; of particular interest is the 
long time behavior of their solutions. In this note we discuss the long time behavior of the solu- 
tions of the stochastic differential equation describing the time evolution of a free quantum particle 
subject to spontaneous collapses in space. We explain why the problem is subtle and report on a 
recent rigorous result, which asserts that any initial state converges almost surely to a Gaussian 
state having a fixed spread both in position and momentum. 



Hilbert space valued stochastic differential equations ap- 
pear all over in quantum physics and their meaning ranges 
from fundamental to effective descriptions of quantum sys- 
tems. They can on the one hand be seen as basic equations 
in collapse models, where the aim is to find a unified de- 
scription of microscopic quantum phenomena and macro- 
scopic classical ones; this is achieved by modifying the 
Schrodinger equation, adding stochastic nonlinear terms 
which model the spontaneous collapse of the wave func- 
tion [1-10]. This fundamental meaning of Hilbert space 
diffusions have been originated in a discrete version, the 
so called GRW model [1], which relies on jump processes 
rather than a continuous diffusion process. One can view 
the diffusion process as continuum limit of the GRW jump 
process [2] ; whether they can be empirically distinguished 
is unclear. Mathematically, diffusion processes are analyt- 
ically easier to handle than jump processes: that is the 
reason why we consider diffusions. 



The same type of equations occur in the theory of con- 
tinuous measurement [11-15]; in this context, they de- 
scribe the effect of continuous measurements on the evo- 
lution of quantum systems. Synonymous to measurement 
is decohcrencc, and therefore the very same effective equa- 
tions appear also in decoherence theory [16,17]. 



The common class of stochastic differential equations 



used in Quantum Mechanics has the following structure: 



dipt 



Hdt + ^f\Y,{Ln-{L n )t)dWi 



-J2(Ln-{L n ) t fdt 



(1) 



with (L n ) t = (ipt\L n ipt) ■ The equation is defined on a suit- 
able domain of a complex separable Hilbert space Tt; it is 
manifestly non-linear, but preserves the norm of the state 
vector. H is the standard quantum Hamiltonian, {L n } is 
a set of commuting sclf-adjoint operators, A is a positive 
constant and {W t } is a family of independent standard 
Wiener process defined on a probability space (f2,JF,P). 
Eq. (1) can be generalized in several different ways, e.g. 
by considering non self-adjoint collapsing operators, nev- 
ertheless the basic structure is preserved. 

One of the most relevant problems in the study of this 
type of equations is to analyze the behavior of their so- 
lutions over characteristic time regimes. The idea, which 
is also been put forward in the context of the classical 
limit of quantum mechanics [18], is that there are three 
time regimes, more or less sharply separated, which char- 
acterize the time evolution of the solution. 1. Collapse 
regime: That is the regime, during which a superposition 
of eigenstates of {L n } collapses to one eigenstate, with a 
probability close to the Born probability rule. This fea- 
ture is a consequence of the martingale structure [4] of 
the non- Schrodinger terms of Eq. (1). In collapse models, 
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the operators {L n } are functions of the position opera- 
tors [2,4,6,8,9], so the wave function collapses in space; 
the constant A is chosen in such a way that, for macro- 
scopic systems, the collapse occurs almost instantaneously. 
In decoherence models instead the collapse time is given 
by the parameters characterizing the interaction with the 
environment [19,20]. 2. Classical regime: The collapsed 
wave packet moves classically; the dynamical equations 
are still stochastic, but the fluctuations around the aver- 
age values — which evolve classically — are very small and 
can be neglected for all practical purposes [2,8]. 3. Dif- 
fusive regime: Eventually, however, the diffusive effect of 
randomness takes over and the third time regime begins; 
the collapsed wave packet starts to diffuse, departing from 
the classical motion. In models of environment-induced 
collapse the diffusion is paired with friction and the classi- 
cal trajectory goes over into a Ornstcin-Uhlenbcck type of 
diffusion. At the end of this note we will come back again 
on these regimes, and will give numerical estimates show- 
ing how they are characterized and separated one from 
the other. The long time regime, which we are concerned 
with in this paper, corresponds to the diffusive one. It is 
the one which can be most easily singled out and treated 
rigorously. 

These time regimes have been studied in particular in 
connection with the following equation, defined in the 
Hilbert space H = C 2 (R): 



dipt 



iih dt + ^ {q ~ {q)t)dWt 

A 



(q-{q)tYdt 



(2) 



q is the position operator and p the momentum opera- 
tor. The equation describes the time evolution of the wave 
function of a free quantum particle subject to spontaneous 
localizations in space. Its importance lies in the fact that 
it is simple enough to be analyzed in full mathematical 
detail; at the same time it is physically interesting and 
gives deep insight in understanding the behavior of more 
complex physical situations. We mention here in partic- 
ular the results of [21] and [10]: in the first case, it has 
been show in mathematical detail how collapse models en- 
sure definite outcomes in measurements on quantum sys- 
tems; in the second case, it has been proven how the entire 
quantum formalism (measurable quantities as self-adjoint 
operators on Hilbert spaces) follows from the basic struc- 
ture of the dynamics of collapse models. In this way the 
dynamical reduction program, initiated with the original 
paper of GRW, finds its completion as a consistent alter- 
native to standard quantum mechanics, at least at the non 
relativistic level. 

One of the most treated problems in connection with 
Eq. (2), as we said, is to establish rigorously the long 
time behavior of its solution [6,8, 12, 16,22-24]. We re- 
port here on a mathematical treatment of this equation 
where the asymptotic behavior is fully and rigorously es- 



tablished [25], correcting gaps in earlier works. At the 
end we will also discuss the three above mentioned time 
regimes in physical units, and we will see how they depend 
on the parameters A and m defining the model. The main 
result of [25] is contained in the following theorem: 

Theorem: Let ip t be the solution of Eq. (2), for the given 
initial condition ip e £ 2 (K), \\tp\\ = 1- Tncn , wit h P- 
probability 1: 



lim ^-Cfe tCt =0, 

t— > + 00 Xt,K t ,C t 



(x - x t ) + ik t x + c t 



with: 

z 2 := (1-i) 
where x t , k t are real stochastic processes, such that 



(3) 

(4) 
(5) 



X + -Kt + V\- fw s ds + X L^W t 
m m l n V m 



-Q{e 



-u,t/2s 



k t = K + V\W t + 0{e- ut/2 ) 



(6) 
(7) 



where X and K are two time-independent random vari- 
ables, while c t is a complex stochastic process and u) = 
2y / h\/m. Of course, c t is such that ^ c is always 
correctly normalized. 

The above theorem states a remarkable property of 
Eq. (2): with P-probability 1 any solution of equation (2) 
converges in the large time limit to a Gaussian solution 
having a fixed spread both in position and momentum, 
while the mean in position x t as well as the mean in mo- 
mentum hk t undergo a random motion guided by W t . 

Instead of elaborating on the mathematical technicali- 
ties which enter the result [25], we would like to explain 
why mathematical rigor is needed to get to a trustable un- 
derstanding. For that, we explain first what goes wrong 
with "intuitive" arguments which have been used in the 
literature, which sound convincing but which nevertheless 
are problematic. 

An argument which has been used in [5-8, 16] to set the 
asymptotic behavior of Eq. (2) and of more general equa- 
tions of the type (1) is based on the large time behavior 
of the variance AA 2 = \\(A — (A) t )i/>t|| 2 of a suitably 
chosen linear operator A, not necessarily self-adjoint. Let 
us assume that the solution ip t of Eq. (1) is such that the 
variance, computed on this state, asymptotically vanishes: 



lim AA 2 = 0. 

t^-\-oo 



(8) 



Since if) is an eigenstate of A if and only if the variance 
AA 2 is zero for such a state, it seems reasonable to con- 
clude that, when (8) holds, ipt converges to an eigenstate 
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of A, which would be the desired result. E.g., in the case 
of Eq. (2) the operator A such that (8) is satisfied is [8,16]: 



A 



2Xm' 



(9) 



which has precisely the coherent states tp^b c defined 
in (4), with a, b £ M, c £ C, as eigenstates. 

The above argument is indeed true when H has finite 
dimension and A is a self-adjoint operator. Let us in fact 
denote by a n its eigenvalues and by (f> n the corresponding 
eigenvectors; for simplicity we assume no degeneracy in 
the spectrum of A. From Eq. (8) one trivially has: 



lim \c n {t)\ 2 [a n ~ (A) t } 2 = 

t— >-+oo 



Vn 



(10) 



withc„(i) = (<f) n \ipt}- Since the coefficients |c„(i)| 2 sum up 
to 1 for any time t, and since (A) t has a unique limit (when 
it exists), then Eq. (10) implies that there exists a n such 
that |c„(t)| 2 — > 6n, n and (A) t — ► <%. As a consequence (j>n 
is the eigenstate towards which ip t converges. 

The above argument fails to apply in the infinite di- 
mensional case. As a counterexample, let us consider any 
orthonormal basis {<fi n : n £ N} of an infinite dimensional 
Hilbert space H, and let A be a self-adjoint operator hav- 
ing <p n has eigenstates relative to the real eigenvalues a n . 
Let us consider the following sequence of vectors: 



+ Put 



+ \Pn\ 2 = l- (11) 



It is easy to show that AA^ = \a n \ 2 \f3 n \ 2 (a n +i — a n ) 2 , 
which asymptotically goes to as soon as (a n +i — o>n) ~ * 
for t — > +oo, without requiring that either a n or fi n 
asymptotically vanishes. In other words, the variance of 
A goes to zero while the tpn need not approach any of the 
eigenstates. 

The reason why the above proof breaks down in infinite 
dimensional spaces is that the spectrum of an operator 
may have an accumulation point, or may be continuous — 
which is actually the case with the operator A defined 
in (9) — in which case the proof for the finite-dimensional 
case does not hold anymore. The physical reason instead 
is that, for such a kind of an operator, the variance, though 
vanishingly small, is always bigger than the difference be- 
tween two eigenvalues, and thus it cannot discriminate 
between the corresponding two eigenstates. 

There is a second argument which has been used to 
study the solution of Eq. (2), in particular its long time 
behavior, which relies on its expansion in terms of coher- 
ent states. It is known [8,13,22-24,26] that Gaussian wave 
functions are solutions of Eq. (2). To see this, let us con- 
sider the following linear stochastic differential equation: 



d<j> t (x) = 



- l T^-dt + V\qdi t -^q 2 dt 
n 2m 2 



( 12 ) 



where £ t is a standard Wiener process defined on a prob- 
ability space (fl,!F,Q). The square norm ||0t|| 2 is a mar- 
tingale [26,27] solving the stochastic differential equation 



d||0 t || 2 = 2^/X(q) t \\4't\\ 2 d^t, and it can be used as a Radon- 
Nikodym derivative of a new probability measure. Let us 
choose Q in such a way that dP/dQ = |j0t|| 2 ; Girsonov's 
theorem implies that the two Wiener processes W t and 
& are related as follows [26,27]: dW t = d£ t - 2V\(q) t dt. 
Given these results, a straightforward application of ltd 
calculus shows that if <j) t solves Eq. (12), than tp t — 
4>t / \\4>t\\ solves Eq. (2). In other words, one can work with 
Eq. (12) — which is linear — and, at the end, normalize the 
wave function and replace £ t with W t according to the pre- 
vious formula in order to get the corresponding solution 
of Eq. (2), for any suitable initial state. 

Direct substitution shows that Gaussian wave functions 



t (x) = exp [-a t (x - a;™) 2 + ik™x + j t ] 



(13) 



solve Eq. (12), if xf, fc™, jt solve suitable stochastic differ- 
ential equations which we do not report here (see e.g. [8,26] 
and references therein), while a t solves the deterministic 
equation: 

d , 2ih 



dt 



A 



m 



(14) 



This equation has the nice feature that a t — > z 2 /2 for 
t — > +oo, for any initial condition a®. In other words, 
any initial Gaussian wave function converges, in the long 
time limit, to a wave function of the type (4). Since C 2 - 
functions can be expressed in terms of coherent states as 
follows: 

<t>(x) = J dx^dk^f(xQ, fcg ) exp [-a (x - x^) 2 + ik„x] , 

(15) 

then using linearity of Eq. (12) and the asymptotic behav- 
ior of a t , one is tempted to conclude that also the general 
solution of Eq. (2) converges to a Gaussian state of the 
type (4). This argument has been used e.g. in [11,23]. 
This conclusion however is false in general. 

That it is false may be seen by scrutinizing the phys- 
ically interesting example of the non-linear but norm- 
preserving differential equation in £ 2 (K): 



Jt Mx) = 



(16) 



The Schrodinger term spreads out the wave function along 
the real axis while the other two terms localize it in space 
(the effect of (q 2 )t is to keep the norm of ipt constant in 
time) ; these two effect are opposite to each other and thus 
cancel when an equilibrium is reached. It is then rea- 
sonable to expect that any initial wave function reaches 
asymptotically such an equilibrium, i.e. it approaches a 
wave function whose spread is the result of the compro- 
mise between the two competing effects. Note that this is 
precisely the kind of intuition concerning the large time 
behavior of the solutions of Eq. (2). 

It is easy to show that a Gaussian wave functions of the 
type (13) are solutions of Eq. (16) if a t solves Eq. (14), 
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while x™, kf and 7 t satisfy the following equations: 



d h A 

— a;? 1 = — A 
ai m 

^ . m a t m 

dt 1 af * 



R * ' 



(17) 
(18) 



the symbols af and a\ denote respectively the real and 
imaginary parts of a t . The two equations for xf and kf 
imply that both these quantities vanish for large times. 
We then have the following result: any initial Gaussian so- 
lution asymptotically converges towards a Gaussian wave 
function centered in the origin of both the position and the 
momentum space, and with a fixed spread equal to z 2 /2. 
We are then precisely in the situation described before, 
which seems to indicate that any solution of Eq. (16) — 
at least those which can be initially written as in (15) — 
converges asymptotically to a Gaussian state with a fixed 
spread. However, we show that this is not true by explic- 
itly finding the general solution of Eq. (16). 

Such an equation can be easily solved by first consider- 
ing the following linear equation: 



ih\q 

2m 



<Mz); 



(19) 



it is immediate to show that if <pt solves Eq. (19), then ipt = 
4>t/\\4>t\\ (assuming that 4> t ^ 0) solves Eq. (16). The above 
equation describes the so called non-self-adjoint (NSA) 
harmonic oscillator and has been already considered in 
the literature [28] , and its general solution is known. The 
eigenvalues of the operator H = p 2 /2m — ihXq 2 are: 



1 - i 









(n-\ 













and the corresponding eigenvectors are: 



u n (x) = yfzexp 



H n (zx), (21) 



where H n {x) is the normalized Hermitc polynomial of de- 
gree n. The general solution of Eq. (19) can then be writ- 
ten as a superposition of the states: 



b?\x) 



exp 



-U) n t 



u n (x), 



(22) 



and the corresponding solution of Eq. (16) can be ob- 
tained by dividing 4>t{x) by its norm; in particular, the 
normalized states ip ( t n) (x) — (f> ( t n) (x)/\\(f> ( t n) \\ are stationary 
solutions of Eq. (16) for any value of n, as the modu- 
lus |$" {x)\ is constant in time: this in particular means 
that they never approach a Gaussian state, with the only 
exception of the ground state ip { t 0} (x) which is already a 
Gaussian function. 

One can easily understand that if <p t is a superposition 
of eigenstates (j> { t n) with n e Af C N, then the normal- 
ized state ipt approaches for large times the stationary 



state -0t"' sucn that n = min{n,n e Af}. In fact, the 
non-normalized eigenstate (j)^ is the one with the weak- 
est time dependent damping factor; accordingly, when one 
normalizes </> t to obtain ip t , this is the only term of the 
superposition for which the damping factor cancels out, 
while all other terms keep an exponential factor propor- 
tional to exp[— (u) n — uOn)t/2\. This means that Eq. (16) has 
not only one stationary state, but infinitely many, since for 
any single-eigenstate solution the damping factor is can- 
celed by the normalization. 

The above analysis disproves the physical intuition 
about the equation and the argument according to which 
properties of Gaussian states are sufficient to analyze the 
behavior of the general solution. The reason why, in this 
case, Gaussian wave functions do not provide a complete 
description is that they always have a non-zero component 
with respect to the ground state ip^ so they necessarily 
converge towards it. But this excludes the situation in 
which the initial wave function does have a null compo- 
nent on such a state. Note also that the states have 
a bigger spread, the bigger the value of n; accordingly, it 
is not even true that Eq. (16) localizes wave functions in 
space, below a fixed spread. 

From a more mathematical point of view, the reason 
why the analysis of the long time behavior in terms of co- 
herent states leads to the wrong conclusion is rather sub- 
tle; we exemplify this by considering the first eigenstate 
Ui(x). Such a state can be expressed as in (15): the Gaus- 
sian term exp [— ao(x — x^) 2 + ik^x\ evolves to a state of 
the form (13), where a tl xf and k$ solve Eqs. (14), (17) 
and (18), respectively, while j t solves the following equa- 
tion: 



d y 
dt 1 * = 



Oif 

1 - 2-r 



(23) 



The function /(xgjfcg 1 ) is highly undetermined, as the 
Gaussian states form an over-complete set; nevertheless it 
has the symmetry property: f{x^ , fcg 1 ) = — /(— x$ , — fcg 1 ), 
since u\(x) = —u\(—x). Then the state at time t reads: 



1} (x) - J dx^dko f(a%, fco ) exp [ 7t ] 

• exp [— a t (x — x™) 2 + ik™x\ ; (24) 



since xf and kf tend to asymptotically, it seems reason- 
able that the second exponential can be easily taken out of 
the integral, so to speak. What is left however is the inte- 
gral of an odd function — as one can see from Eq. (23), the 
global factor 7 t is even in the initial conditions x^, fcg 1 — 
which of course is zero. Therefore the large time behavior 
requires much more scrutiny and cannot be read off from 
the "zeroth-order approximation" . 

The above discussion is useful not only because it shows 
that care must be taken when analyzing the large time 
behavior of the solutions of Eq. (2), but also because it 
helps in understanding what actually happens in this case. 
As proven in [25] each solution solution of Eq. (12), from 
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which the corresponding solution of Eq. (2) can be directly 
obtained, can be mapped into a solution of Eq. (19) and 
can be written as follows: 



<t>t{x) = e 



+ OC 

n=0 



■'"UC»> 



(x-x t ), (25) 



where 4> t (x) are exactly the eigenstates of the NSA har- 
monic oscillator defined in (22), while x t and k t are the 
stochastic processes introduced in (6) and (7), and j t is 
a new stochastic process function of c t appearing in (4); 
moreover 



+ OC 

X] ak + r 

k=0 



V(k + n)l 



(26) 



£ t being another stochastic process; the complex terms a m 
are the coefficients of the superposition of the initial state 
(f>o in terms of the eigenstates </> m) = u m . 

A crucial point is that the series defined in (25) is norm 
convergent only for t > (4c + l)/w, where c is a constant. 
This means that the representation of the general solu- 
tion of Eq. (12) given in (25) is suitable only for studying 
the large time behavior of tpt solution of Eq. (2), not its 
properties at finite small (compared to u>) times, which are 
also very interesting — even more interesting than the long 
time behavior, from a physical point of view — as discussed 
in [25]. 

As we can see from (25), the collapse mechanism is ba- 
sically the same as in the previous example of Eq. (16): 
ipt = 0t/||0t|| converges to the lowest non- vanishing eigen- 
state appearing in the sum (25), which is the only one 
whose exponential damping factor is fully canceled by nor- 
malization. However, a crucial difference with respect to 
the solutions of Eq. (16) occurs, which is embodied in the 
form (26) of the coefficients a 1 ^ : whichever the initial 
condition, due to the Brownian motion the solution <p t of 
Eq. (12) always picks a non-vanishing component on the 
ground state (i.e. a t ^ a.s., for any t > 0); accord- 
ingly, the normalized solution ip t converges almost surely 
towards such a state. This is how the asymptotic conver- 
gence of the general solution of Eq. (2) towards a Gaussian 
state occurs. The detailed analysis can be found in [25]. 

Discussion of the time regimes. We conclude this 
note by showing how one can read from (25) the three time 
regimes we mentioned at the beginning, which depict the 
main features of the evolution of the solution i[) t of Eq. (2). 
As discussed e.g. in [8], it is physically convenient to take 
A proportional to the mass m of the particle, according to 
the formula: 

A = A — , (27) 
too 

where too is a reference mass and Ao is the universal cou- 
pling constant which applies to all systems. Accordingly, 
to is the only free parameter of the model. The above as- 
sumption embodies the fact that, for a composite system, 



the coupling constant Acm associated to the motion of the 
center of mass of the systems turns out to be the sum of 
the coupling constants Aj associated to the motion of each 
constituent. In models of wave function collapse — which 
are the only ones we will refer to, in the following - in or- 
der for the collapse mechanism to have the same strength 
of that of the original GRW model [1], one has to set 
Ao ~ 1.00 x 10~ 2 m -2 sec -1 when too is taken equal to 
the mass of a nucleon (~ 1.67 x 10~ 27 Kg). 

The first important point to stress is that the series 
in (25) is convergent only for t > 1 := (4c + l)/w; since c 
is of order 1 and u ~ 5.01 x 10~ 5 sec -1 (in a GRW-like 
scenario), then t ~ 2.00 x 10 4 sec. As we shall soon see, 
this implies that the representation (25) is valid only for 
the classical and diffusive regime, to which the long time 
behavior belongs. If one wants to analyze the collapse 
regime and possibly the onset of classical motion other 
representations of the solution must be used, e.g. the one 
in terms of the Green's function. Given this, here we shall 
use (25) to describe all three time regimes, keeping in mind 
that, for t <t, the representation (25) makes sense only 
if a finite sum of terms appears. This is nevertheless suf- 
ficient to describe, qualitatively, the time evolution of the 
solution of Eq. (2). 

According to Eq. (22), the n-th term of the series (25) 
defining <f> t decays exponentially in a time equal to l/nu>, 
which means that after a given time T all terms with 
n > 1/ivT have decayed. This is the collapse regime, 
previously introduced: any initially spread out wave func- 
tion collapses in space. However, according to (20) the 
frequency uj depends on the ration A/to and thus is inde- 
pendent of the mass of the particle: this means that both 
for micro- and for macro-objects the eigenstates decay in 
the same way; in particular, in a GRW-like model, after a 
time T ~ 10~ 3 sec, which is about the perception time of 
a human being, only the states with n > 2.00 x 10 7 have 
decayed. This is a very large number, which seems not to 
justify the reason why, according the collapse models, the 
wave function collapses the faster, the bigger the system. 

However, the parameter z 2 defined in (4), which sets the 
width of the eigenstates u n of the NSA harmonic oscillator 
through Eq. (21), grows linearly with to. This implies that 
the spread in position a™ of the state u n decreases with 
the mass according to the formula: 



(28) 




nass = mn 



Accordingly, after a given time T, the states with n < 
1 /ujT which have survived the collapse, all have a spread, 
which is the smaller, the bigger the system, and thus sum 
up to form a wave function whose extension, for fixed T, 
decreases as to increases. This is how the amplification 
mechanism works, which says that the collapse regime is 
shorter, the bigger the mass of the system. This behavior 
can be seen also by inspecting directly Eq. (2): if one 
performs the replacement x — > y = y / mjrn x, then the 
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dependence of the equation on m disappears. This means 
that the only effect of the mass is to "shrink" the wave 
function in space, not to influence directly its collapse in 
time. The fact that the collapse occurs in good agreement 
with the Born probability rule is not so easy to see; ref. [25] 
contains a discussion of this point. 

Let us now consider the case of a macroscopic particle 
and let us suppose that the wave function has collapsed to 
a function which, on macroscopic scales, is almost point- 
like. Then, within macroscopic scales, the average value 
(q)t and (p) t of the position and momentum operators, 
respectively, are practically equal to xt and hk t defined 
in (6) and (7), i.e. with the average position and momen- 
tum of the Gaussian ground state. We can then analyze 
x t and hk t in place of (q) t and (p) t ; e.g. let us consider 
the time evolution of x t . 

As we see from Eq. (6), the deterministic part of the evo- 
lution corresponds to the solution of Newton's laws, while 
the random part of the evolution is proportional to ra -1 / 2 , 
thus it is suppressed for large s ystem s. E.g.: VXh/m ~ 
2.57 x 1CT 19 m sec -3 / 2 and jh~f^{ ~ 3.24 x 1(T 16 m 
sec -1 / 2 , for a 1-g object. The figures show that for very 
long times the (almost point-like) wave function moves 
practically along a straight line, like a classical particle. 
This corresponds to the classical regime, where the mo- 
tion is essentially classical. This regime last the longer, 
the bigger the system; for macroscopic systems, as the nu- 
merical example shows, it is much longer than the time 
during which the system can be kept isolated, so that its 
dynamics can be described by Eq. (2). 

The classical regime ends when the fluctuations become 
relevant on a macroscopic scale, i.e. essentially when Wt ~ 
\Jm/h which is about 9.53 x 10 10 sec if we take Angstrom 
units for length. After this time, the fluctuations start 
to blur the deterministic motion: this is what we have 
called the diffusive regime, which dominates the remaining 
evolution of the wave function. 
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